Dynamical interferences to probe short-pulse photoassociation of Rb atoms 

and stabilization of Rb 2 dimers 
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We analyze the formation of Rb 2 molecules with short photoassociation pulses applied to a cold 
85 Rb sample. A pump laser pulse couples a continuum level of the ground electronic state X X E^" with 
bound levels in the 0^ (5S+5P!/ 2 ) and 0+ (5S+5P 3 / 2 ) vibrational series. The nonadiabatic coupling 
between the two excited channels induces time-dependent beatings in the populations. We propose 
to take advantage of these oscillations to design further laser pulses that probe the photoassociation 
process via photoionization or that optimize the stabilization in deep levels of the ground state. 

PACS numbers: 32.80.Qk, 33.80.Ps, 34.50.Rk 



Making ultracold molecules in the lowest vibrational 
level v = of the ground electronic state and creating 
stable molecular condensates is presently an important 
challenge since it opens the road toward ultracold chem- 
istry l2fl • Schemes based on photoassociation (PA) of 
ultracold atoms [|[ with cw lasers, have been very suc- 
cessful to form molecules in an excited electronic state. 
The latter have been stabilized into excited vibrational 
levels of the ground electronic state fi@, but not yet into 
v = except for the case of RbCs p . The possibility of 
controlling PA by use of short laser pulses has been dis- 
cussed in theoretical papers @, H, H, El and very recently 
attempted by two experimental groups [ll|, [Tl] , both in 
the rubidium case. Success in such experiments will cre- 
ate a bridge between the two domains of cold matter and 
coherent control, where femtosecond (fs) pulses are used 
to control chemical reactions [l3|. Unfortunately, up to 
now, PA experiments with fs laser pulses have achieved 
destruction of the molecules already existing in the trap 
rather than creation of additional molecules [ll], [12] ■ 

Finding ways to avoid this destructive effect is there- 
fore a crucial step in the development of experiments. 
A promising route is PA through the resonant cou- 
pling mechanism as realized, with cw lasers, in Cs2 [H[ 
and RbCs (la |. For the case of Rb2, it populates the 
0+(5S+5P 1/2 ) and 0+(5S+5P 3/2 ) coupled series. This 
is a textbook example of global mixing of two molecu- 
lar vibrational series due to spin-orbit (SO) coupling and 
manifested by strong perturbations in the Rb 2 0+Jluores- 
cence and cw photoassociation spectra (cf. Refs. |mll7|). 
The aim of the present paper is to draw attention to the 
coherent character of the time evolution in the coupled 
excited states after the PA pulse, and to the possibility or 
taking advantage of the subsequent dynamical interfer- 
ences to optimize stabilization and avoid the destruction 
of the molecules. We consider a simple model with pulses 
in the picosecond range which populate a limited number 
of bound vibrational levels, and we analyze characteristic 
time-dependent oscillations that should appear in exper- 
imental signals. 

Time-dependent fluorescence signals manifesting the 
coupling between deeply bound levels of the two se- 
ries have been previously observed by pump-probe spec- 



troscopy [la ]. In this experiment, a molecular beam of 
Rb 2 was excited by a 120-fs laser pulse from the v = 
level of X : E+ =X to low (~ 5600 cm" 1 below the 5S+5P 
asymptote) vibrational levels of the A : E+ =A state, with 
a classical inner turning point ~ 8oq. A delayed pulse, 
operating around the outer turning point of the b 3 II u =b 
potential (12ao), probes the population transferred to 
this state. It presents an oscillatory behavior with sev- 
eral characteristic periods of the order of 1 ps. This is 
explained in terms of interferences between several paths 
for the motion of the wavepacket in the excited state due 
to the crossing, at i? s hort ~ 9.3ao, of the potential energy 
curves (hereafter referred to as Va and V&), correspond- 
ing to the two Hund's case a A and b states. 

This crossing can be seen in Fig. [Ha) , where we plot 
Va and V& as determined in Ref. [19j . We also show 
the ground-state potential obtained by matching ab initio 
calculations 120(1 to the long-range dispersion potential 
— J2 n C n /R n [21(, and adjusting the repulsive barrier to 
reproduce the large value of the scattering length [22j . 

In this paper, we consider PA at a detuning of 
« 60 cm" 1 below the Di line (1 cm" 1 w 30 GHz), 
exciting loosely bound vibrational levels of the coupled 
system in the vicinity of the up-to-now unexplored cross- 
ing at long range i?i on g ~ 29.2ao [Fig. Q^c)]. The Va 
curve crosses the V& potential curve shifted down from 
Vb by A-Ef s /3 = 79 cm -1 by the inclusion of the diagonal 
part of the spin-orbit effective Hamiltonian, Hso [1 91 ] • 

The diagonalization of Hso within the {A, b} sub- 
space, hereafter referred to as diabatic, renders two 
Hund's case c 0„ (Pi/ 2 ) and 0+(P 3 / 2 ) adiabatic curves, 
correlated with the 5S+5P!/ 2 and 5S+5P 3 / 2 asymptotes, 
split by AE[ S . The mixing angle 9(R) defines the trans- 
formation from the diabatic to the adiabatic basis at a 
given internuclear distance R: 

|0+(P V a)> = cos9(R) |A) + sm9(R) |b) , (la) 
|0+(P 3 /a)> =-sia9(R) \A)+cos9(R) |b> . (lb) 

In the adiabatic representation, the two excited channels 
are coupled by radial coupling, governed by the deriva- 
tive d9/dR. We display in Fig. QJb) the i?-dependence 
of cos 2 9. There is a sharp variation, corresponding to a 
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FIG. 1: (color online) (a) Potential energy curves for the 
ground state X 1 E^" (dash-dotted line), Va (black, solid) and 
Vt, [red (gray), solid], (b) |A) component in |0u(Pi/2)) as 
given by Eqs. fT}. (c) Zoom of the energy range where PA 
takes place, showing both diabatic {A, b} (solid lines) and 
adiabatic (dashed) potentials. The shaded box indicates the 
PA window in the energy domain for the pulse studied. The 
dotted lines are at i? s h or t and -Ri ong , and symbols stand for 
the classical turning points of the levels contributing to the 
decomposition of the excited wavepacket in the uncoupled 0^ 
basis. 



singlet-triplet change of character, in the region of -R s hort, 
where the two diabatic curves cross abruptly. In the 
range of distances from 10an to 15an, the splitting be- 
tween the diabatic curves Va and Vj, is very large com- 
pared to AE{s, and 0+(Pi/ 2 ) has a pure singlet character. 
For R > 15a , cos 2 9{R) decreases toward its asymptotic 
value 1/3, which, due to the very similar slopes of Va 
and Vb at large R [cf. Fig. Hfc)], is reached only at dis- 
tances far beyond 40arj. In contrast with the crossing at 
-Rshort explored in Ref. [18|, this long-range crossing is 
not localized. 

We perform calculations for a system of 85 Rb atoms at 
a temperature of 100 fiK as in the usual photoassociation 
experiments [H|,[l2]- The 85 Rb isotope is chosen since the 
manifestation of the resonant coupling between the two 
series is more remarkable [17] . Only s-wave scattering 
and J = rotational levels are considered for simplicity, 
even though a full treatment of the rotational structure 
will be required to compare with experimental data when 
they become available. We consider a chirped pulse of 
duration tq , that delivers an energy -E pu isc uniformly over 
an area a. It is centered at time tp and has a frequency 
that varies linearly in time, u>(t) = W£+X-(i— tp), around 
the carrier frequency lol. The laser is red detuned from 



the atomic Di line at wdi by = H{u>t>i — w^); x is the 
linear chirp rate in the time domain. The instantaneous 
intensity I(t) of the pulse involves a Gaussian envelope 
with a full width at half maximum equal to tq [H|. For 
this pulse, 98% of the energy E pu \ se is delivered in the 
time window [tp—Tc, tp + Tc] [8j, the instantaneous laser 
frequency being then resonant with all the excited levels 
with a binding energy in the range [Al — h\\\Tq, Al + 
^lxl r c]i which defines a PA window in energy [231 ] . cf. 
Fig. [He). 

We choose a nonperturbative description for the dy- 
namics of the PA process and of the vibration of 
the molecules, and solve the time-dependent coupled 
Schrodinger equations [1, 0] to compute the wavepacket 
motion in the ground X and the excited A and b 
states. The X and A channels are coupled by M7(£) = 
— y/I(t)/(2ceo)V(R), where c and £n are the speed of 
light and permittivity of vacuum, and ~D(R) is the molec- 
ular transition dipole moment; at large distances D(R) = 
4.245ean (e = proton charge, an = Bohr radius). A ra- 
dial grid up to 30 000an has been built with the mapped 
grid method [24| to faithfully represent the initial state 
as a stationary scattering level in X, and the vibrational 
wavefunctions in the excited states. The time propaga- 
tion was performed as in Ref. Q. 

From now on, we concentrate on a short PA pulse, red- 
detuned by w 60 cm -1 from the Di line (corresponding 



to A 



pump 



798 nm). The bandwidth is set to 15 cm 



large enough to resonantly couple the initial state |X, E = 
98.85 /xK) with two levels of the 0+(P 3 / 2 ) series and 13 
levels of the 0+(Pi/ 2 ) series, but small enough to avoid 
population of continuum levels in the excited state [cf. 
Fig. HJc)]. The pulse is centered at tp — 50 ps, with 
t c = 10 ps, x = 4.41 x 10~ 2 ps~ 2 , and E pu i sc = 41 nJ 
focused on a = 2.8 x 10~ 3 cm 2 . We choose x > to 
maximize the population transfer (25|. 

We present in Fig. [2] the evolution of the popula- 
tion in the excited states, \^ exc (R,t)\ 2 = \^ A (R,t)\ 2 + 
^(i?,*)! 2 , well after the PA pulse has finished. The PA 
probability for a single pair of atoms is 4.21 x 10~ 6 , and 
most excited population concentrates around two peaks 
at 22an and 37arj. The relative importance of these peaks 
changes in time. For example, the amplitude of the peak 
at 22arj oscillates with a period T* « 8 ps. 

This beating results from the interferences between 
the population of the stationary levels v' of the cou- 
pled system of excited states. Indeed, in the present 
case, the coupling is such that neither the diabatic 
nor the adiabatic stationary vibrational levels consti- 
tute a good basis for the analysis. Here, the coupled 
states |</>™ up ), which have components on the two elec- 
tronic channels, are labeled by v' according to their in- 
creasing energy E v i, In this basis, the decomposition 
of the two-component wavepacket reads \I' cxc (i?, t) — 



,coup 



(R) exp(—iE v it/H), where each c v > as well as 



the total excited population are constant after the pulse. 

We show in Figs. G3h) and (3Jc) the radial density 
|</>™ up (i?)| 2 of two stationary states present in the ex- 
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FIG. 2: (color online) Evolution of the wavepacket in the 
excited channels ^ cyic (R,t) after the PA pulse has finished. 
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FIG. 3: (color online) Stationary wavefunctions of (a) the 
Rb^~(X 2 £+) level bound by 264 cm" 1 ; (b) the 0+ level 
54 cm -1 below Di; and (c) the 0^ level 51 cm -1 below Di. 



cited wavepacket. The lower panel corresponds to a level 
51 cm -1 below the 5S+5Pi/ 2 limit. It turns out to be as- 
cribable to a "0+ (Pi/ 2 ) level", as most of its probability 
density is accumulated close to the outer turning point 
-Rext,i/2 ~ 37an, corresponding to this potential and this 
energy [black circle in Fig. QJc)]. The state in Fig.[3jb), 
which corresponds to a level 54 cm" 1 below Di , is a typ- 
ical example of a "resonant 0+(P 3 / 2 ) level" [l7|: it has a 
55% weight in the adiabatic P3/2 channel, and features 
two maxima of probability, corresponding to the turning 
points i? cx t,i/2 and -R ox t,3/2 ~ 22a [red (gray) circle in 
Fig. Die)] in 0+(P 1/2 ) and 0+(P 3/2 ) respectively. As the 
instantaneous frequency of the chirped pulse is resonant 
with ~15 levels in the coupled 0+ series, ^ exc will have 
components on stationary levels of both types, which ex- 
plains the double-peak structure observed in Fig. 

Now, it is easy to show that the population density at 
an internuclear distance R has a time dependence deter- 
mined by all the beating frequencies cjy = (J3j — Ej)/ H, 
where E^ and Ej are the energies of those levels for 
which cv 7^ 0. The largest weights in the decomposi- 
tion correspond to the levels of the coupled basis labeled 
by v' = 406 and 408, two "0+(P 1/2 )" levels separated 
by AE ss 4.24 cm" 1 , which corresponds to a beating of 
Tbcat ~ 7.85 ps, in agreement with T* observed in Fig. [2] 

The observed oscillation in the population density close 
to -R e xt,3/2 can be probed experimentally with a suit- 
able laser pulse that ionizes the molecule and coherently 
populates several levels of the lowest Rb^X 2 !]^) po- 
tential [26j with classical turning point i? cx t,+ ~ 22an. 
The stationary wavefunction of one such level is plot- 
ted in Fig. [3]Ja). Assuming that the ionization process 
is a vertical transition, the population transfer to the 
ionic channel will concern the part of the wavepacket 
that is close to i? ex t,-H i-e., the peak around i? ox t,3/2- 
Indeed, we show in Fig. [Ha) the time evolution of the 
overlap £„„ |(Rb+(X 2 S+), </'|*exc(*))| 2 for the levels 
v" e [185, 194] that have i? xt,+ ~ (21 - 23)a and 
can be populated using a laser with central wavelength 
Aprobc = 479.8 nm and bandwidth w 30 cm" 1 . (The over- 
lap is calculated with both A and b components of the 
excited wavepacket, as no selection rules apply to the ion- 



izing transition.) The signal has a rich structure of peaks 
and troughs. In the Fourier transform (FT) of a longer- 
duration signal [cf. Fig. Etc)], one observes a characteris- 
tic time T\ = 6.2 ps (and its multiples), which is identi- 
fied with the beating period between the levels v' = 407 
[a "0+(P 1/2 )" level] and v' = 410 [a "resonant 0+(P 3/2 )" 
level, cf. Fig. [3^b)]. Thus, T\ is a characteristic time of 
the P1/2-P3/2 population transfer. Quite remarkably, a 
different ionizing pulse designed to probe the wavepacket 
close to 37an (A pro be = 469.3 nm) would render a sig- 
nal dominated by T 2 = 2Tb ca t = 15 ps, the beating pe- 
riod between the two "0+(P 1/2 )" levels v' = 407,408 [cf. 
Fig. 0|b) and0Jc)]: as the classical turning point of the 
0+ (P3/2) components is 22an, where the RbJ vibrational 
wavefunctions have no appreciable amplitude, almost no 
effect of the 0+(Pi/ 2 )-0„ (P 3 / 2 ) coupling is expected to 
show up at this distance. Of course, some peaks are visi- 
ble in Fig. |3Jc) that correspond to other beating periods 
cj,^ 1 . Nevertheless, T\ and T 2 can be used as fingerprints 
of the two different ionization processes, as each of them 
appears only in one signal. 
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FIG. 4: (color online) Evolution of the overlap of |>]/exc(t)) 
with (a) 10 Rb+(X 2 £+) levels with i? cxt ,+ » 22a ; (b) 16 
Rb^(X 2 E^") levels with -R ex t,+ ~ 37ao. (c) Fourier transform 
of signals in (a) [blue (black) line] and (b) [red (gray) line] 
in the time window t £ [70, 300] ps. The dashed vertical line 
stands for tp — 50 ps and arrows indicate the most relevant 
timescales. 
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FIG. 5: (color online) Evolution of the overlap of |$ e xc(f)) 
with (a) 17 Rb 2 (X 1 E+) levels with E hia d < 24 cm" 1 ; (b) 28 
Rb 2 (X 1 E+) levels with E hind w 1300 - 2500 cm -1 . 



We next analyze the dynamics after a dump pulse 
toward bound vibrational levels of Rb2(X 1 E+). Note 
that due to selection rules, this transition is selecting 
the A 1 !]^ component in the excited wavepacket. The re- 
sults are shown in Fig. [H^a) for the vibrational levels with 
binding energy £bind ^ 24 cm" 1 , and in Fig. [5jb) for the 
levels with 1300 cm -1 < £bind < 2500 cm" 1 , which have 
a reasonable Franck-Condon factor with $ cxc and are 
reachable by a fs dump pulse (Adump = 708 nm). In both 
cases, Tbeat and the mean value of the vibrational pe- 
riod in the pure 0+ (P1/2) potential in this energy range, 
T-2 ~ 13 ps, have an important role, and can be used to 
determine the best timing for the dump pulse. 

In summary, we have studied the dynamics of Rb2 0+ 
molecules created by a picosecond PA pulse from cold 
85 Rb atoms. The excited wavepacket spans ~ 15 vibra- 
tional levels in the coupled 0+ basis and presents two 
maxima of probability at i? e xt,i/2 and i? C xt.3/2- The sub- 
sequent dynamics shows quantum interferences, with a 
beating of the population close to -R C xt.3/2- This beat- 



ing can be monitored by photoionizing the wavepacket 
with a laser with well-defined energy spectrum, which in 
practice defines a window of internuclear distances whose 
density probability is probed. The corresponding signal 
and its Fourier transform can serve as identification tool 
for ongoing PA experiments [TTI . [I2I ] . We have also ana- 
lyzed a pump-dump pair of pulses to form deeply-bound 
ground-state molecules. The resulting time dependence 
would allow, in an experiment with a sequence of pump- 
dump pairs at a high repetition rate, to choose the delay 
between consecutive pairs, so that a pump pulse would 
not dissociate the molecules created by the preceding 
dump pulse. Moreover, these molecules could be fur- 
ther transferred to even deeper levels by other pairs of 
pulses (see also Ref. [HJ). We studied here a pulse with 
a relatively large detuning to avoid coupling to the con- 
tinuum levels in the excited channels. In the future, for 
a closer comparison with experiments, femtosecond PA 
pulses with smaller dctunings and a mask to cut the blue 
part of the pulse spectrum that transfers populations to 
the continuum, will be analyzed. Also, we will further 
study the role of SO coupling in 87 Rb and 133 Cs sam- 
ples. 
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